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ABSTRACT 

We present radial velocity measurements for 74 globular clusters (GCs) in the nearby 
giant elliptical NGC 5128, of which 31 are newly discovered clusters. All the GC can¬ 
didates were taken from the list of possible new clusters given in the Harris, Harris, 
& Geisler (2004) photometric survey. In addition to the newly confirmed clusters, we 
identified 24 definite foreground stars and 31 probable background galaxies. From a 
combined list of 299 known GCs in NGC 5128 with measured radial velocities and 
metallicity-sensitive ( C — T±) photometric indices, we construct a new metallicity dis¬ 
tribution function (MDF) for the cluster system. The MDF shows an approximately 
bimodal form, with centroids at [Fe/H] = —1.46 and —0.53, and with nearly equal num¬ 
bers of metal-poor and metal-rich clusters in the two modes. However, there are many 
intermediate-color clusters in the distribution, and the fainter clusters tend to have a 
higher proportion of red clusters. These features of the MDF may indicate a widespread 
age range within the cluster system as well as an intrinsically broad metallicity spread. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 
5128) — globular clusters: general — techniques: radial velocities 


Wisiting Observer, Cerro Tololo Interamerican Observatory, Operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. 
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1. Introduction 

Globular clusters (GCs) are simple stellar populations (single age, single metallicity) that help 
to trace the formation history of their host galaxy. GCs are also powerful tools in revealing the 
halo kinematics and dynamics, if large samples of globular cluster radial velocities can be obtained. 

NGC 5128, the central galaxy in the Centaurus cluster 3.9 Mpc away, and the nearest readily 
accessible giant E galaxy, provides a unique and exciting possibility to study the dynamics and for¬ 
mation history through its globular cluster system (GCS). However, finding its clusters individually 
is a difficult job: at its moderately low Galactic latitude of 19°, field contamination by both fore¬ 
ground stars and background galaxies is significant, and the GCs cannot be distinguished by only 
image morphology, object size, or color. Thus except for the very expensive and time-consuming 
approach of resolving the GCs into individual stars by (say) HST imaging, radial velocity determi¬ 
nations provide the only definitive method for classifying a candidate object in the NGC 5128 field 
as a genuine GC. 

NGC 5128 has a statistically determined total population of 980 ±120 globular clusters, giving 
it a low specific frequency S/v = 1.4±0.2 (Harris, Harris, & Geisler 2004, hereafter HHG04). An Sn 
value this low (in the range occupied by large spiral galaxies) is suggestive of a merger origin (Harris 
2001); other indicators of mergers are the faint isophotal shells in the halo of the galaxy (e.g. Peng 
et al. 2002), and the activity in the inner 5 kpc region indicated by the visible dust lane, the recent 
star formation, and the presence of gas. The GCS, with its wide mixture of cluster metallicities 
and perhaps ages, should provide many additional clues to the formation history surrounding this 
giant elliptical. An excellent starting point for any such discussion is to have as large as possible a 
list of individual, certainly identified GCs through radial velocity measurements. 

Previous work has led to a list of 215 confirmed GCs in NGC 5128 with radial velocity mea¬ 
surements (Peng, Ford, & Freeman 2004, hereafter PFF04). Many hundreds more clusters remain 
to be found. In this paper, the first of a new series of studies of the GCS in this unique galaxy, we 
present new velocity measurements and the discovery of 31 new GCs. 


2. Data Reduction and Velocity Measurements 

The wide-held photometry study of the NGC 5128 held in the Washington CMT\ system by 
Harris et al. (2004) provides a list of 327 GC candidates having globular-cluster-like colors and 
slightly nonstellar morphology. This list provides the basis for our new radial velocity survey. 
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2.1. Observations and Data Reduction 

Observations were taken at the Cerro Tololo Inter-American Observatory (CTIO) 4m telescope 
in April 2004. The Hydra multi-object fibre-fed spectrograph, with 40' diameter held of view, was 
used with the KPGL1 grating and 2x binning in dispersion, giving a spectral resolution element of 
approximately 3.5 binned pixels (FWHM). The spectral coverage was A3900 A — A5300 A, extending 
from Ca II H+K at the blue end to the Mgb features at the red end at 1.18 A/px. Of our three 
assigned nights, two were lost to clouds, and the third had poor seeing throughout {2" — 3"). 
Nevertheless, we successfully observed two Hydra held configurations, both centered on the galaxy. 
Total exposure times were 6 x 1800 sec for the hrst held and 7 x 1800 sec for the second. 

The data reduction was performed through the dohydra package within IRAF 1 , following the 
normal sequence of preprocessing described on the Hydra webpages. The raw images were overscan 
and bias subtracted, and the detector was hat-helded with milk hats (daytime sky exposures) using 
a diffusing hlter in front of the Hydra hbres. Cosmic rays were subtracted with a Laplacian edge 
detection method (van Dokkum 2001). Spectral throughput and dispersion corrections were applied 
based on the standard continuum lamp and Hydra HeNeAr penray lamp to obtain wavelength 
calibrated spectra. 


2.2. Background Light Removal 

Removal of background light from each object spectrum includes both the uniform dark-sky 
illumination and (for the GC candidates close to the galaxy center) the bulge light from NGC 5128 
itself. To trace both effects, we placed 40-50 hbers in each Hydra held conhguration on blank-sky 
positions sprinkled across the entire held, then plotted the mean intensity of each sky hber versus 
radius from galaxy center. The data showed that the bulge light profile was well matched by a 
(1/r 2 ) radial dependence (see, e.g. Harris, Harris, & Geisler 2004). 

For radii larger than r ~ IT, the galaxy bulge light contamination was negligible and we simply 
subtracted the average of the sky spectra at these larger radii from all the object spectra. For the 
inner objects, we then removed the additional bulge light contamination through appropriate radial 
scaling of the average bulge sky spectrum. 


2.3. Radial Velocities 

From the fully preprocessed and background-subtracted spectra, we derived radial velocities 
obtained through Fourier cross-correlation, using the fxcor task in the radial velocity IRAF package. 


1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy Inc., under cooperative agreement with the National Science Foundation. 
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The spectrum of each candidate object was correlated against a template spectrum of known 
velocity and the correlation peak was fit with a Lorentzian function. Our template spectrum was 
constructed from the sum of five of the previously known bright GCs in NGC 5128 (C7, Cll, C17, 
C19, and C20), which we observed in both of our Hydra fields and which have well determined 
velocities taken from PFF04. Thus, the template was very well tuned to the the spectral features 
in the GCs that we are searching for since it was built from exactly the same type of object. 


3. Results 

3.1. Candidate Distribution 

NGC 5128 has a systemic velocity of 541 km s -1 (Hui et al. 1995) and an internal stellar 
velocity dispersion of ~ 140 km s _1 (Wilkinson et al. 1986). These features allow a clean velocity 
separation between GCs and foreground stars in the Milky Way or background galaxies. As shown 
(e.g.) by PFF04, the foreground stars have velocities almost all less than ~ 200 km s -1 , while any 
faint background galaxies will have velocities of many thousands of km s _1 . Any objects that fall 
in the range 200 < v r < 1000 km s” 1 (and drawn from the candidate list pre-selected by image 
morphology, magnitude and color) will therefore be near-certain GCs. 

From both Hydra fields combined, we obtained radial velocity measurements for 86 candidate 
GCs and 43 previously known GCs. The previously known clusters were used to help establish the 
zero point of our velocity scale and check the internal measurement uncertainties (see below). Of 
the candidates, we found 31 with 240 < v r < 900 km s _1 , which we identify as newly discovered 
GCs. Of these 31, 5 are X-ray sources previously detected by Kraft et al. (2001). In Table 1, 
we list the data for these GCs. Successive columns give the identification number (WHH prefixes 
denote our new objects), coordinates, radial distance from the center of NGC 5128, V and (C — T\) 
photometry from HHG04, radial velocity, and velocity uncertainty from the fxcor solution. 

Of the 55 remaining candidates, 24 are definite foreground stars with —91 < v r < 110 kms -1 , 
and 31 are likely background galaxies; these are listed in Table 2. The objects with null results for 
the radial velocity are ones whose spectra had a poor correlation with the GC template spectrum, 
suggestive of a background galaxy with either very large redshift or spectral features very different 
from the template. Confirming 31 candidates as GCs out of the 86 candidates observed yields an 
encouraging 36% hit rate given the extremely high field contamination level, and demonstrates the 
value of rigorous preselection by color and object morphology (see the discussion of HHG04). 

Table 3 lists the resulting radial velocity measurements of the 43 previously known GCs re¬ 
measured in this study. The identification number prefixes C, G, or pfLgc are originally listed by 
Harris et al. (1992) and PFF04. 

Figure 1 shows the locations of the newly confirmed and previously known GCs. Whereas most 
previous surveys have concentrated on fields along the isophotal major axis of NGC 5128 (along 
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a line running SW to NE), it is worth noting that many of our newly identified GCs lie along the 
less well studied minor axis parallel to the inner dust lane. There are likely to be many more to 
be found still. Figure 2 shows the velocity distribution of the GCs and foreground stars measured 
from this study; the two types of objects are clearly distinguished by radial velocity at v r ~ 200 
km s -1 . 

The analysis of the GCS radial profile by HHG04 showed that the majority of GCs lie within 
r = 15' from the galaxy center, with the main contribution of the clusters residing in an even more 
compact environment within r = 10b Our new sample of 31 GCs has a median value of r = 8.25b 

The mean radial velocity difference between the previously known clusters in our study and 
those of PFF04 is less than 10 km s _1 , well within the standard deviation of the mean (see below). 
The zero point of our velocity scale therefore needs no adjustment. 


3.2. Comments on Individual Objects 

The raw spectra of clusters WHH-12, WHH-25, WHH-31, pfLgc-073, and the non-GC candi¬ 
dates 100, 133, and 136 had a few intensity spikes that could not be removed in the preprocessing 
stage. These regions were artificially removed before the cross-correlation solutions, resulting in 
higher than normal velocity uncertainties (see Tables 1 & 2). 

For the previously known cluster CIO, we measure v r = 855 ± 33 km s^ 1 . This differs quite 
a bit from the PFF04 measurement of 522 ± 53 km s -1 , but it agrees well with Hesser, Harris, & 
Harris (1986) (770 km s _1 ) and also with the new measurement from Beasley et al. (2005) (836 km 
s^ 1 ). We therefore believe our measurement to be in the correct range. 

There is one candidate for which the velocity solution revealed a second correlation peak only 
slightly smaller than the highest peak that is automatically found and fit by the fxcor task. For 
WHH-22, the highest correlation peak yields v r = 492 ± 126 km s -1 , clearly classifying this object 
as a GC; but a second peak yields v r ~ 20 km s -1 , which would classify it as a foreground star. 
Visual inspection of the candidate spectrum shows hydrogen and calcium line strengths consistent 
with those for clusters. We keep this object as a probably GC, but with these uncertainties in 
mind. A strong secondary peak was also revealed for the known cluster, pff_gc-060. The highest 
peak gives v r = 818 ± 96 km s -1 , but a second peak gives 217 km s -1 , which would put it on the 
borderline between GCs and foreground stars. Its spectrum does not have high enough signal-to- 
noise to determine upon visual inspection whether or not it is a cluster. However, this object has 
been measured by PFF04 to have v r = 893 ± 19 km s -1 , agreeing with our highest peak correlation 
value. 

For three previously known GCs, C21, pfLgc-089, and pfLgc-064, the highest correlation peak 
fit generated by fxcor yielded radial velocities that were inconsistent with previous measurements. 
A closer look at the correlated spectra shows a second strong correlation peak in all cases. Choosing 
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the secondary peak gave v r = 461 ±96 km s _1 for C21 (compared with 465 ±7 km s _1 from PFF04) 
and v r = 553 ± 103 km s _1 for pff_gc-089 (compared with v r = 502 ± 61 km s _1 from PFF04). 
C21 has also been previously measured by Hesser, Harris, & Harris (1986) to be at v r = 490 km 
s^ 1 . For these two objects we therefore adopted the second correlation peaks. The known cluster, 
pff_gc-064, measured in both field configurations, presents a slightly different problem: the initial 
poorly correlated fit in field 1 nominally classified it as a possible background galaxy. Yet, there is 
a second strong correlation peak with v r = 548 ± 97 km s -1 , agreeing with the field 2 result for the 
same object of 563 ± 84 km s -1 and also with the PFF04 measurement of v r = 557 ± 31 km s -1 . 

There are three other candidates measured in both field configurations, 43/44, 109/110, and 
101. Their radial velocities classified these objects as possible background galaxies, with their 
averaged radial velocities as displayed in Tables 2. Objects 43 and 44 in the candidate list of 
HHG04 turn out to be the same object, and similarly candidates 109 and 110 are the same. Lastly, 
we note that 12 of the HHG04 candidate objects measured in this study turned out to be previously 
known GCs. These objects are listed in Table 4. 


3.3. Velocity Measurement Uncertainties 

The radial velocity uncertainties we quote in the tables are from the fxcor solutions. Generally, 
these are well correlated with the object’s apparent magnitude and thus the signal-to-noise of the 
spectra. Figure 4 shows the nominal radial velocity uncertainties plotted versus V magnitude. The 
new GCs (excluding WHH-12, WHH-25, and WHH-31; see above) with V < 18.5 had a mean 
(j{v r ) = ±32 km s -1 and S/N = 13 per pixel in the continuum region near A4525A For V > 18.5, 
the mean uncertainty was a(v r ) = ±74 km s^ 1 and S/N (A4525A) = 4.5 per pixel. 

We tested the fxcor uncertainties in several ways. First, the comparison between our velocities 
and those of PFF04 for 42 clusters in common (excluding CIO as mentioned above) is shown in 
Figure 3. The standard deviation of the differences between the two studies was ±44.5 km s -1 . If 
the two studies have comparable internal uncertainties, then this suggests that for GCs brighter 
than V — 19 the typical uncertainty of our measurements is no larger than ±30 km s _1 , in agreement 
with the mknoise simulations in Fig. 4 (see below). 

Another internal test was to compare the velocities for objects that we measured in both Hydra 
field settings. There were 13 such GCs measured twice, 12 of them brighter than V = 19: C2, 
C6, C13, C14, C18, C19, C20, C26, C37, C38, C44, pfLgc-056, and pff_gc-064 (V = 19.80 mag). 
For these, the mean velocity difference between the two fields was 3.9 km s -1 , with a standard 
deviation of ±22.3 km s -1 , again consistent with Fig. 4. 

Finally, we compared the fxcor uncertainties with those generated from simulated noisy spec¬ 
tra, adding Gaussian noise to the template spectrum with the IRAF artdata.mknoise task. Several 
dozen of these artificial spectra were cross-correlated against the original template spectrum and 
the resulting radial velocity uncertainties plotted against the S/N for each noisy spectrum. The 
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results verify that for S/N > 5 per pixel (corresponding roughly to V < 18.5), the velocity uncer¬ 
tainty was smaller than ±30 km s^ 1 . Fainter than this, the uncertainty rises exponentially, reaching 
±100 km s -1 at S/N ~ 2 per pixel (or V 20). In Fig. 4, the solid line shows the velocity uncer¬ 
tainty increasing with V magnitude, from the mknoise simulations. This ideal line is the minimum 
uncertainty due only to random noise, providing a lower boundary to the observed uncertainties. 

In summary, we find that the results generated by fxcor (with its uncertainties as listed in 
Tables 1, 2, & 3) provide estimates of the internal velocity uncertainties that are quite close to the 
estimates from simulations and from comparisons with previous data. This may in large part be 
because we deliberately used a standard velocity template which accurately matches the program 
objects. We conclude that the internal uncertainties of our measured velocities are near ±20 km 
s^ 1 for the very brightest clusters, increasing to ±100 km s _1 for the ones at V ~ 20 or fainter. 


4. The Metallicity Distribution Function 

A new analysis of the radial velocities and kinematics of the confirmed GCS in NGC 5128, 
with a much enlarged cluster sample over that used in PFF04, will be discussed in a later paper 
(Woodley et al. 2005). For the present, we provide here only a brief discussion on the metallicity 
distribution function (MDF). 

The top panel in Figure 5 shows the MDF for a total of 299 GCs in NGC 5128 that are 
radial-velocity members and that have available (C — Tj) photometric indices: 168 are previously 
confirmed clusters (PFF04), 31 are the new clusters from this study, and 100 more are from the 
upcoming study of Beasley et al. (2005). Although there are a few dozen other highly probable 
GCs that have recently been found in NGC 5128 from image morphology studies (Rejkuba 2001; 
Harris et al. 2002; Gomez et al. 2005), we restrict the current discussion to an almost certainly 
“pure” sample generated by the combination of object color, morphology, and radial velocity. 

The dereddened colors {C — Ti)o from HHG04 were transformed to [Fe/H] through the conver¬ 
sion derived by Harris & Harris (2002) calibrated through the Milky Way cluster data. We adopt a 
foreground reddening value of E(B — V) = 0.11 for NGC 5128, corresponding to E(C — T\) = 0.22. 
Because the transformation is slightly nonlinear, the uncertainty in [Fe/H] is a function of metal¬ 
licity. For a typical color uncertainty a(C — Tf) ~ 0.1, the corresponding uncertainty in [Fe/H] is 
near ±0.2 dex at the metal-poor end, decreasing to ±0.07 dex at the metal-rich end. 

The MDF, binned in steps of A [Fe/H] = 0.10, is shown in Figure 5, along with the same 
distribution for the Milky Way clusters (data from Harris 1996). As has been found in all previous 
studies of NGC 5128 starting with Harris et al. (1992), the MDF is extremely broad and populated 
by roughly equal numbers of metal-poor and metal-rich clusters. The color bimodality of NGC 5128 
has recently been confirmed by Peng, Ford, & Freeman (2004) in 4 different photometric indices. 
They used spectroscopic line indices of H/3 versus [MgFe]' to suggest the blue GC population was 
nretal-poor and old (8 — 15 Gyr) and the red GC population was metal-rich and young (mean age 



~ 5+2 Gyr). A bimodal-Gaussian distribution (shown in Fig. 5) provides a reasonable fit to the 
total MDF, but not as cleanly as in many other giant galaxies including the Milky Way (plotted 
for comparison in the lower panel of the figure). The fit of the bimodal model to the MDF is not 
well constrained, allowing a wide range of particular combinations of Gaussians. The two centroids 
shown, at [Fe/H] = —1.46 and —0.53, are similar to the fit obtained by HHG04 (see their Figure 10) 
from a smaller cluster sample. The metal-poor part of our entire sample makes up ~ 54% of the 
total population and the metal-rich part ~ 46%, in close agreement with HHG04. It does not seem 
likely that internal differential reddening of the clusters is an important factor in creating the color 
spread, since reddening-generated differences in color of several tenths of a magnitude would be 
needed to move a cluster from the metal-poor mode to (e.g.) [Fe/H] ~ — 1, yet all the objects with 
(C — T\ ) photometry lie clearly outside the dust lane and into the clear regions of the bulge and 
halo (see Harris et al. 2004, for extensive discussion). 

As has long been realized, the Milky Way MDF (bottom panel in Figure 5) provides a clean 
match to the bimodal-Gaussian model, with peaks at [Fe/H] = —1.58 ± 0.05 and —0.64 ± 0.07. 
The standard deviations are 0.32 and 0.23 for the metal-poor and metal-rich subpopulations, both 
narrower than their counterparts for NGC 5128. Both centroids are slightly but noticeably more 
nretal-poor than those in NGC 5128. 

By comparison, the MDF of the M31 globular clusters shows [Fe/H] peaks that are more 
similar to NGC 5128. Perrett et al. (2002) find M31 centroids of [Fe/H] = —1.44 ± 0.03 and 
—0.50 ± 0.04. Thus, it is likely that if NGC 5128 is a product of major mergers, then it would 
need M31-sized galaxies as progenitors, rather than systems like the Milky Way which have lower- 
nretallicity clusters. 

If we assume that NGC 5128 was built through a series of mergers extending for several 
Gigayears, then the GCS might contain a wide range of ages. Thus the MDF shown here, based on 
photometric metallicities from (C — T\) with the assumption that all the clusters are homogeneously 
old, would underestimate the true metallicity of younger clusters. This effect would artificially 
increase the spread in the MDF and partially blur out a distribution that was intrinsically bimodal. 
The classic age/metallicity degeneracy for broad band colors prevents us from distinguishing this 
alternative from an intrinsically broad, even MDF. 

Figure 6 displays the NGC 5128 sample divided into bright (V < 19.5) and faint (V > 19.5) 
clusters. Interestingly, the fainter sample appears to have a higher proportion of metal-rich clusters 
than the brighter sample. A direct least-squares correlation of [Fe/H] versus V magnitude shows 
no significant trend within the larger scatter in [Fe/H], It is possible that we are simply seeing 
an intrinsically broad and uniform spread of cluster metallicities, in which the younger and more 
metal-rich ones tend to be of somewhat lower mass or luminosity, consistent with the age differences 
suggested by PFF04. A wide range of ages would further increase the apparent spread of the MDF 
as seen in Fig. 6. 

Further spectroscopic studies and more detailed modelling will be needed to test these spec- 
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illations. In later analyses (Woodley et al. 2005; Beasley et al. 2005) we will discuss the age 
distributions and kinematics of the cluster system further. 

In Figure 7, the cluster metallicities are displayed as a function of projected galactocentric 
radius. Both metal-poor ([Fe/H] > —1.0) and metal-rich ([Fe/H] < —1.0) clusters have been found 
extending out to a radius of 40' (~ 45 kpc). However, beyond r ~ 20', there are double the 
number of metal-poor compared to metal-rich clusters. The mean metallicity beyond this distance 
is [Fe/H] = —1.21 with a standard deviation of 0.47. As is discussed in more detail in HHG04, this 
change in ratios of subpopulations is not likely to be due to selection bias and is therefore probably 
intrinsic to the system. 


5. Summary 

We present new radial velocities for 86 globular cluster candidates selected from Harris, Harris, 
& Geisler (2004). The results show 31 newly identified clusters with 240 < v r < 900 kms 1 , 24 
definite foreground stars with —91 < v r < 110 km s _1 , and 31 probable background galaxies. We 
also obtained velocities for 43 previously known GCs in NGC 5128. The internal uncertainty in 
our v r measurements ranges from ±20 km s -1 at the bright end of the sample (V < 18) to ~ ±100 
km s -1 at V ~ 20. 

Combining our data with previously published material, we derive a metallicity distribution 
for the GCs based on 299 clusters. A bimodal distribution with centroids at [Fe/H] = —1.46 and 
—0.53 provides an approximate fit to the MDF, though not highly accurate, and the data hint that 
the relative proportions of metal-poor and metal-rich clusters may change with luminosity. Two 
possible interpretations are that (a) we may be seeing the combined effects of a several-Gy internal 
age spread in the cluster system, combined with a trend for the younger clusters to be progressively 
less massive, or (b) the MDF is genuinely broad and evenly distributed, unlike the more cleanly 
bimodal MDFs seen in other giant galaxies. 


We greatly appreciate the data reduction advice and support provided by Knut Olsen at 
CTIO, and also advice regarding uncertainty tests from David James at Vanderbilt University. We 
thank Mike Beasley for transmitting cluster velocity data in advance of publication. This work was 
supported by the Natural Sciences and Engineering Research Council of Canada through research 
grants to GLHH and WEH. 
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Table 1. Radial Velocity Measurements for Globular Clusters 


Cluster 

RA 

(J2000) 

Dec 

(J2000) 

Rgc 

(arcmin) 

V 

(mag) 

(C-Tr) 

V r 

(km s' 1 ) 

a(v r ) 

(km s' 1 ) 

WHH-1 

13 

24 

21.40 

-43 

02 

36.8 

12.19 

18.23 

1.484 

600 

35 

WHH-2 

13 

24 

23.98 

-42 

54 

10.7 

13.56 

19.87 

1.270 

582 

81 

WHH-3 

13 

24 

32.17 

-43 

10 

56.9 

14.10 

19.59 

1.470 

636 

78 

WHH-4 

13 

24 

40.60 

-43 

13 

18.1 

14.89 

19.13 

1.927 

685 

43 

WHH-5 

13 

24 

44.58 

-43 

02 

47.3 

8.04 

19.50 

1.386 

671 

41 

WHH-6 

13 

24 

47.37 

-42 

57 

51.2 

8.06 

19.75 

1.931 

712 

67 

WHH-7 a 

13 

25 

05.02 

-42 

57 

15.0 

5.68 

17.43 

1.863 

722 

20 

WHH-8 a 

13 

25 

07.62 

-43 

01 

15.2 

3.66 

18.05 

2.026 

690 

32 

WHH-9 

13 

25 

08.51 

-43 

02 

57.4 

3.93 

18.90 

1.979 

315 

100 

WHH-10 

13 

25 

12.84 

-42 

56 

59.8 

4.95 

19.70 

1.334 

664 

141 

WHH-1 l a 

13 

25 

14.24 

-43 

07 

23.5 

6.71 

19.55 

2.066 

562 

55 

WHH-12 

13 

25 

18.27 

-42 

53 

04.8 

8.25 

19.31 

1.521 

558 

99 b 

WHH-13 

13 

25 

21.29 

-42 

49 

17.7 

11.91 

19.34 

1.854 

467 

47 

WHH-14 

13 

25 

25.49 

-42 

56 

31.2 

4.64 

20.28 

1.300 

436 

80 

WHH-15 

13 

25 

26.78 

-42 

52 

39.9 

8.48 

20.13 

1.153 

506 

57 

WHH-16 a 

13 

25 

27.97 

-43 

04 

02.2 

2.89 

19.20 

2.036 

749 

80 

WHH-17 

13 

25 

29.25 

-42 

57 

47.1 

3.38 

18.73 

1.409 

588 

127 

WHH-18 

13 

25 

30.07 

-42 

56 

46.9 

4.39 

18.48 

1.752 

771 

31 

WHH-19 

13 

25 

31.03 

-42 

50 

14.9 

10.92 

17.51 

1.177 

470 

49 

WHH-20 

13 

25 

34.36 

-42 

51 

05.9 

10.12 

19.19 

1.517 

242 

54 

WHH-21 

13 

25 

35.22 

-43 

12 

01.5 

10.97 

19.22 

0.994 

243 

61 

WHH-22 a 

13 

25 

35.31 

-43 

05 

29.0 

4.56 

18.63 

1.631 

492 c 

126 

WHH-23 

13 

25 

45.90 

-42 

57 

20.2 

5.07 

18.60 

1.137 

286 

63 

WHH-24 

13 

25 

46.00 

-42 

56 

53.0 

5.43 

19.85 

1.288 

566 

48 

WHH-25 

13 

25 

50.34 

-43 

04 

08.2 

5.12 

20.17 

1.621 

551 

95 b 

WHH-26 

13 

25 

56.59 

-42 

51 

46.6 

10.76 

19.16 

1.762 

412 

36 

WHH-27 

13 

26 

12.82 

-43 

09 

09.2 

11.51 

18.73 

1.885 

545 

60 

WHH-28 

13 

26 

14.18 

-43 

08 

30.4 

11.25 

19.32 

1.414 

506 

123 

WHH-29 

13 

26 

22.08 

-43 

09 

10.7 

12.79 

19.98 

1.538 

505 

78 

WHH-30 

13 

26 

23.60 

-43 

03 

43.9 

10.55 

19.56 

1.488 

470 

66 

WHH-31 

13 

26 

41.43 

-43 

11 

25.0 

16.96 

19.65 

1.618 

573 

67 b 
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Table 1—Continued 


Cluster RA 

(J2000) 

Dec 

(J2000) 

Rgc 

(arcrnin) 

V (C-Ti) 

(mag) 

v r 

(km s -1 ) 

(j(v r ) 

(km s” 1 ) 



a Detected as X-ray source from Kraft et al. (2001) 

b Regions of spectrum affected by cosmic ray spikes were masked out during correla¬ 
tion. 

c Strong secondary correlation peaks. 







Table 2. Foreground Stars and Background Galaxies 


ID 

R.A. 

Dec. 

V r 

cr(ry) 

Comments 


(J2000) 

(J2000) 

(km s” 1 ) 

(km s _1 ) 



36 

13 

24 

02.55 

-42 

49 

29.4 

64 

42 

star 

3 

13 

24 

28.96 

-43 

10 

58.3 

-24 

23 

star 

8 

13 

24 

52.11 

-43 

08 

24.3 

-1 

26 

star 

19 

13 

24 

56.46 

-43 

06 

34.7 

-46 

24 

star 

11 

13 

25 

04.20 

-43 

02 

29.3 

-30 

26 

star 

14 

13 

25 

08.59 

-43 

07 

43.4 

-17 

47 

star 

13 

13 

25 

10.06 

-43 

07 

40.7 

-30 

21 

star 

26 

13 

25 

17.79 

-43 

13 

05.6 

-70 

42 

star 

25 

13 

25 

24.37 

-43 

12 

23.1 

-13 

23 

star 

62 

13 

25 

25.10 

-43 

12 

27.0 

-82 

79 

star 

27 

13 

25 

42.62 

-43 

10 

21.4 

-20 

32 

star 

77 

13 

25 

47.98 

-43 

00 

57.6 

106 

106 

star 

4 

13 

25 

55.77 

-43 

03 

00.9 

18 

21 

star 

1 

13 

25 

56.40 

-42 

59 

37.3 

40 

23 

star 

2 

13 

25 

57.38 

-42 

59 

48.5 

1 

22 

star 

6 

13 

25 

58.39 

-42 

55 

27.6 

-22 

21 

star 

23 

13 

25 

58.44 

-43 

03 

06.4 

-3 

25 

star 

20 

13 

26 

10.41 

-43 

00 

53.6 

-11 

26 

star 

73 

13 

26 

15.41 

-42 

57 

46.8 

3 

118 

star 

56 

13 

26 

17.81 

-43 

00 

29.7 

-60 

37 

star 

17 

13 

26 

20.74 

-43 

09 

13.2 

51 

33 

star 

89 

13 

26 

21.11 

-43 

03 

08.2 

-85 

99 

star 

61 

13 

26 

32.79 

-43 

02 

59.4 

31 

50 

star 

38 

13 

26 

33.08 

-43 

03 

14.0 

-91 

40 

star 

283 

13 

24 

04.50 

-42 

48 

16.2 

47636 

81 

galaxy 

285 

13 

24 

05.49 

-42 

47 

24.1 

4729 

124 

galaxy 

290 

13 

24 

13.26 

-42 

57 

42.1 

- 

- 

PC C 

94 

13 

24 

28.15 

-42 

53 

04.6 

35382 

178 

galaxy 

67 

13 

24 

28.46 

-43 

14 

56.3 

- 

- 

PC 

122 

13 

24 

29.70 

-43 

02 

06.4 

- 

- 

PC 

;/44 

13 

24 

30.86 

-42 

52 

40.8 

35527 a 

70 

galaxy 






Table 2—Continued 


ID 

R.A. 

Dec. 

v r 

a(v r ) 

Comments 


(J2000) 

(J2000) 

(km s -1 ) 

(km s -1 ) 



117 

13 

24 

34.38 

-42 

57 

32.4 

- 

- 

PC 

127 

13 

24 

34.61 

-43 

12 

50.5 

7717 

194 

galaxy 

133 

13 

24 

35.30 

-43 

06 

15.4 

- 

_b 

PC 

119 

13 

24 

48.71 

-42 

52 

35.5 

- 

- 

PC 

298 

13 

25 

01.09 

-42 

48 

21.9 

- 

- 

PC 

136 

13 

25 

01.16 

-43 

11 

59.0 

- 

b 

PC 

299 

13 

25 

01.29 

-42 

54 

26.5 

2544 

92 

galaxy 

132 

13 

25 

23.45 

-42 

53 

26.2 

2017 

78 

galaxy 

108 

13 

25 

25.73 

-43 

05 

16.6 

- 

- 

PC 

103 

13 

25 

28.72 

-42 

50 

12.3 

- 

- 

PC 

314 

13 

25 

37.53 

-42 

55 

33.5 

- 

- 

PC 

116 

13 

25 

39.63 

-43 

04 

01.4 

32124 

253 

galaxy 

123 

13 

25 

41.37 

-43 

12 

12.6 

31018 

230 

galaxy 

317 

13 

25 

45.63 

-43 

01 

15.5 

5850 

91 

galaxy 

72 

13 

25 

49.23 

-43 

00 

02.2 

7349 

87 

galaxy 

90 

13 

25 

49.26 

-43 

02 

20.4 

9942 

113 

galaxy 

135 

13 

25 

49.95 

-43 

02 

26.3 

38298 

114 

galaxy 

107 

13 

25 

56.63 

-42 

58 

45.5 

- 

- 

PC 

101 

13 

26 

02.24 

-43 

08 

55.7 

_ a 

- 

PC 

109/110 

13 

26 

02.23 

-43 

17 

34.8 

_ a 

- 

PC 

134 

13 

26 

17.27 

-43 

06 

39.3 

24742 

79 

galaxy 

113 

13 

26 

21.29 

-42 

57 

19.1 

- 

- 

PC 

125 

13 

26 

23.64 

-43 

00 

45.6 

- 

- 

PC 

100 

13 

26 

50.44 

-42 

50 

36.4 

- 

_b 

PC 


a Measured in both fields; the radial velocity is the average. 

b Regions of spectrum ignored during correlation to remove spikes, leading 
to high uncertainties. 

C PC denotes a poor correlation of the object to the GC template spectrum. 
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Table 3. Radial Velocity Measurements for Previously Known Globular Clusters 


Cluster 

RA 

Dec. 

Rgc 

V (C-Ti) 

v r 

a(v r ) 


(J2000) 

(J2000) 

(arcmin) 

(mag) 

(km s -1 ) 

(km s _1 ) 


C2 

13 

24 

51.47 

-43 

12 

11.2 

12.86 

18.38 

1.546 

606 a 

38 

C3 

13 

24 

58.21 

-42 

56 

10.0 

7.33 

17.68 

1.940 

559 

20 

C4 

13 

25 

01.81 

-43 

09 

25.5 

9.53 

17.88 

1.451 

687 

49 

C6 

13 

25 

22.19 

-43 

02 

45.6 

1.89 

16.83 

- 

846 a 

18 

C7 

13 

26 

05.40 

-42 

56 

32.4 

8.30 

16.90 

1.533 

599 

15 

CIO 

13 

24 

48.04 

-43 

08 

14.3 

10.13 

18.30 

1.727 

856 

33 

Cll 

13 

24 

54.70 

-43 

01 

21.6 

6.02 

17.67 

2.011 

765 

15 

C13 

13 

25 

06.22 

-43 

15 

11.6 

14.58 

18.48 

1.703 

599 a 

25 

C14 

13 

25 

10.49 

-42 

44 

52.8 

16.56 

17.75 

1.655 

691 a 

19 

C15 

13 

25 

30.39 

-43 

11 

49.6 

10.69 

18.43 

1.881 

672 

33 

C17 

13 

25 

39.72 

-42 

55 

59.1 

5.62 

17.49 

1.422 

785 

18 

C18 

13 

25 

39.86 

-43 

05 

01.8 

4.48 

17.29 

1.603 

476 a 

19 

C19 

13 

25 

43.38 

-43 

07 

22.9 

6.87 

17.87 

0.324 

607 a 

18 

C20 

13 

25 

49.68 

-42 

54 

49.3 

7.50 

17.89 

1.472 

745 a 

16 

C21 

13 

25 

52.73 

-43 

05 

46.5 

6.52 

17.62 

1.576 

461 

96 

C26 

13 

26 

15.25 

-42 

48 

29.4 

15.36 

17.94 

2.068 

370 a 

21 

C30 

13 

24 

54.35 

-42 

53 

24.7 

9.84 

17.05 

1.788 

786 

17 

C32 

13 

25 

03.37 

-42 

50 

46.1 

11.29 

18.22 

2.006 

708 

32 

C33 

13 

25 

16.26 

-42 

50 

53.3 

10.47 

18.67 

- 

520 

41 

C36 

13 

26 

07.73 

-42 

52 

00.2 

11.72 

18.24 

1.378 

696 

43 

C37 

13 

26 

10.57 

-42 

53 

42.7 

10.81 

18.17 

1.691 

612 a 

21 

C38 

13 

26 

23.77 

-42 

54 

01.1 

12.50 

18.34 

1.544 

397 a 

23 

C39 

13 

26 

42.00 

-43 

07 

44.9 

15.11 

17.28 

1.647 

247 

26 

C44 

13 

25 

31.73 

-43 

19 

22.8 

18.25 

18.48 

1.441 

515 a 

43 

C47 

13 

25 

49.94 

-42 

52 

09.3 

9.88 

18.80 

1.402 

592 

34 

C48 

13 

25 

49.83 

-42 

50 

15.1 

11.63 

18.65 

1.281 

557 

77 

C50 

13 

26 

19.64 

-43 

03 

18.6 

9.75 

18.59 

1.633 

605 

90 

G271 

13 

25 

13.94 

-42 

57 

42.5 

4.25 

18.80 

1.092 

263 

98 

G342 

13 

25 

05.75 

-42 

59 

00.3 

4.53 

18.20 

- 

565 

37 

G369 

13 

24 

57.49 

-42 

59 

23.3 

5.78 

18.80 

1.436 

527 

58 

;c-018 

13 

24 

47.08 

-43 

06 

01.7 

8.87 

18.91 

1.603 

554 

44 






Table 3—Continued 


Cluster 

RA 

Dec. 

R S c 

V (C-TO 

v r 

a(v r ) 


(J2000) 

(J2000) 

(arcmin) 

(mag) 

(km s -1 ) 

(km s" 1 ) 


pff_gc-041 

13 

25 

11.14 

-43 

03 

09.6 

3.62 

19.61 

1.471 

432 

179 

pff_gc-042 

13 

25 

12.19 

-43 

16 

33.7 

15.67 

18.97 

1.352 

678 

87 

pff_gc-056 

13 

25 

32.80 

-42 

56 

24.3 

4.84 

18.64 

1.286 

275 a 

44 

pff_gc-058 

13 

25 

35.12 

-42 

56 

45.1 

4.60 

18.68 

1.259 

376 

78 

pff_gc-060 

13 

25 

42.42 

-42 

59 

02.5 

3.43 

18.93 

1.634 

818 c 

96 

pff_gc-064 

13 

25 

43.89 

-42 

50 

42.5 

10.85 

19.80 

1.758 

556 a,c 

63 

pff_gc-073 

13 

25 

52.76 

-42 

58 

41.6 

5.21 

20.10 

1.759 

470 

189 b 

pff_gc-075 

13 

25 

53.47 

-43 

03 

56.6 

5.49 

19.61 

1.768 

764 

59 

pff_gc-076 

13 

25 

53.74 

-43 

19 

48.5 

19.26 

19.07 

1.923 

351 

39 

pff_gc-083 

13 

26 

01.81 

-42 

58 

15.0 

6.89 

20.04 

1.683 

473 

53 

pff_gc-085 

13 

26 

06.39 

-43 

00 

38.1 

7.11 

19.63 

1.262 

469 

108 

pff_gc-089 

13 

26 

20.19 

-43 

10 

35.7 

13.48 

18.95 

1.324 

553 

103 


a Measured in both fields; the radial velocity listed is the average. 

b Regions of spectrum affected by cosmic ray spikes were masked out during correlation. 
c Strong secondary correlation peaks. 
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Table 4. Previously Known Clusters in HHG04 Candidate List 


Candidate ID Cluster 


9 

C18 

12 

C3 

16 

C14 

30 

C2 

42 

pfLgc-056 

48 

pfLgc-060 

82 

pfLgc-041 

83 

pfLgc-075 

85 

pff_gc-085 

95 

pfLgc-064 

111 

pff_gc-083 

118 

pff_gc-073 
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Fig. 1.— The positions in RA and Dec: (left) the 31 new globular clusters found in this study; 
(right) 215 known globular clusters from previous literature (Peng, Ford, and Freeman 2004). The 
cross hairs indicate the centre of the galaxy. 
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ae velocity distribution of the 74 measured GCs (v r ~ 200 to 900 
—90 to 110 km s -1 ). 
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Fig. 3.— Comparison of our radial velocity to the PFF04 radial velocity for the 42 known clusters 
in common (excluding CIO; see text). The standard deviation of the differences is ±44.5 km s _1 . 
The average internal uncertainty of each set of measurements is shown in the upper left of the plot. 
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Fig. 4.— The fxcor velocity uncertainties obtained for the 74 globular clusters listed 
& 3, compared to their V magnitudes. The solid line is the ideally expected increasi 
uncertainty with increasing V magnitude, based on the mknoise simulations. 
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Fig. 5.— Top: The metallicity distribution function for 299 GCs including the 31 new GCs from 
this survey. The distribution is fit with two Gaussians with centroids at [Fe/H] = (—1.46, —0.53) 
and standard deviations (0.54, 0.30). Bottom: The metallicity distribution function for 148 Milky 
Way GCs with Gaussian centroids at [Fe/H] = (—1.584=0.05,—0.64T0.07) and standard deviations 
(0.32, 0.23) 
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Fig. 7.— Metallicity from the C-Ti colour index vs. projected galactocentric distance for the 
radial-velocity confirmed GCs. The crosses denote clusters from Harris et al. (1992); van den 
Bergh, Hesser, &; Harris (1981); Hesser, Harris, & Harris (1986), the solid circles are clusters from 
Peng, Ford, & Freeman (2004), and the open circles are new clusters from this study. 




